arXiv:1507.07654vl [cond-mat.str-el] 28 Jul 2015 


DMFT+E approach to disordered Hubbard model 


E. Z. Kuchinskii, M. V. Sadovskii 
Institute for Electrophysics, Russian Academy of Sciences, Ural Branch, 

Amundsen str. 106, Ekaterinburg 620016, Russia 
M.N. Mikheev Institute for Metal Physics, Russian Academy of Sciences, Ural Branch, 

S. Kovalevskaya str. 18, Ekaterinburg 620990, Russia 

We briefly review the generalized dynamical mean-field theory DMFT-I-E treatment of both repul¬ 
sive and attractive disordered Hubbard models. We examine the general problem of metal-insulator 
transition and the phase diagram in repulsive case, as well as BCS-BEC crossover region of attrac¬ 
tive model, demonstrating certain universality of single - electron properties under disordering in 
both models. We also discuss and compare the results for the density of states and dynamic con¬ 
ductivity in both repulsive and attractive case and the generalized Anderson theorem behavior for 
superconducting critical temperature in disordered attractive case. A brief discussion of Ginzburg 
- Landau coefficients behavior under disordering in BCS-BEC crossover region is also presented. 
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I. INTRODUCTION 


Strongly correlated electronic systems, which are 
mainly realized in a range of compounds contain¬ 
ing transition or rare-earth elements with partially 
filled 3d, 4/ or 5/ shells, attract attention of sci¬ 
entists because of their unusual physical properties 
and are notorious for major difficulties in theoret¬ 
ical description. Perhaps the most significant de¬ 
velopment in this area was the discovery of high 
temperature superconductivity in copper oxides, 
which are considered to be the typical example of 
strongly correlated systems. 

Early qualitative ideas formulated mainly by 
Mott [y as well as the introduction of the seminal 
Hubbard model 0 inspired the hundreds of theo¬ 
retical papers, which now constitute the separate 
branch of condensed matter theory. Probably the 
most impressive achievement if this field in recent 
years was the development of dynamical mean-field 
theory (DMFT), which provides an asymptotically 
exact solution for the Hubbard model in the limit 
of infinite dimensions [^-Q. 

Most of the studies of strongly correlated sys¬ 
tems within Hubbard model are devoted to the 
case of repulsive interactions among electrons, 
which are directly related to many topical prob¬ 
lems, with most attention payed to the physics 
of high-Tc superconductivity in cuprates and the 
general problem of metal-insulator transition in 
cuprates and other similar oxides of transition met¬ 
als. 

Another direction of research is the studies of the 
Hubbard model with attractive interaction, which 
is related mainly to rather old problem of strong 
coupling superconductivity, especially to the the¬ 
oretical description of the notorious BCS to BEC 
(Bardeen - Cooper - Schrieffer to Bose - Einstein 
Condensation) crossover, which is also directly re¬ 
lated to the problem of high-Tc superconductivity 
in copper oxides. Starting with pioneering papers 


by Eagles and Leggett at T = 0 and impor¬ 

tant prog ress achieved by Nozieres and Schmitt- 
Rink [n|, who suggested an effective method to 
study the transition temperature crossover region, 
this field has produced the large number of the¬ 
oretical papers published during the recent years, 
including the successful! applications of DMFT ap¬ 
proach. 

This last area of research is also directly con¬ 
nected with recent progress in experimental studies 
of quantum gases in magnetic and optical dipole 
traps, as well as in optical lattices, with control¬ 
lable parameters, such as density and interaction 
strength (cf. reviews [il,[i3|), which has increased 
the interest to superconductivity (superfluidity of 
Fermions) with strong pairing interaction, includ¬ 
ing the region of BCS - BEC crossover. 

In recent years we have develo ped the so called 
generalized DMFT-bS approach |l4l - [l7l| . which is 
very convenient for the studies of different addi¬ 
tional interactions in repulsive Hubbard model, 
such as pseudogap fluctuations [T3 - [l7l| . disorder 
BEl, electron - phonon interaction HI) etc. 
This approach is also well suited to analyze two- 
particle properties, such as optical (dynamic) con¬ 
ductivity [iMllIl- In Ref. 221 we have used this 
approximation to calculate single - particle prop¬ 
erties of the normal phase and optical conductivity 
in attractive Hubbard model. Recently DMFT-I-E 
approach was used by us to study disorder influ¬ 
ence upon superconducting transition temperature 
in this model [1^ . 

Below we shall concentrate on discussion the of 
disorder effects in both repulsive and attractive 
Hubbard models. There are not so many works, 
devoted to the studies of disorder effects in Hub¬ 
bard models, because of many theoretical com¬ 
plications, related to the problem of mutual in¬ 
terplay of disorder scattering and Hubbard inter¬ 
action. We shall concentrate exclusively on our 
DMFT-I-S approach, which is actually very conve¬ 
nient here and provides good interpolation scheme 
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between different limiting cases. We shall discuss 
the results obtained in our previous work, similar¬ 
ities and dissimilarities of disorder effects in both 
repulsive and attractive Hubbard models, demon¬ 
strating in certain cases the universal dependences 
on disorder. 



FIG. 1: Typical “skeleton” self-energy diagramms in 
DMFT-I-S approximation. 


II. THE BASICS OF DMFT+E 
APPROACH IN DISORDERED SYSTEMS 

The Hamiltonian of disordered Hubbard model 
can be written as: 

H = -t'^ + + H ^71^1714, ( 1 ) 

{ij )(7 i <7 i 

where t > 0 is the transfer integral between near¬ 
est sites of the lattice, U is the onsite interaction 
([/ > 0 in the case of repulsive interaction, while in 
the case of attraction U < 0), rua = the 

operator of the number of electrons on the lattice 
site i, Qia (ajcr) i® annihilation (creation) operator 
for electron with spin a on site i. The local energy 
levels Ci are assumed to be independent random 
variables at different lattice sites (Anderson dis¬ 
order) To simplify diagram technique in the 
following we assume the Gaussian distribution of 
these energy levels: 

Parameter A represents here the measure of disor¬ 
der and this Gaussian random field (with “white 
noise” correlation on different lattice sites) gener¬ 
ates “impurity” scattering and leads to the stan¬ 
dard diagram technique for calculation of the en¬ 
semble averaged Green’s functions [ 2 ^. 

Generalized DMFT-I-E abroach [l3 - [T^ ex¬ 
tends the standard DMFT introducing an 

additional self-energy Sp(e) (in general case mo¬ 
mentum dependent), which is due to some inter¬ 
action mechanism outside the DMFT. It gives an 
effective procedure to calculate both single- and 
two-particle properties mull- The single-particle 
Green’s function is then written in the following 
form: 

G'(£,p) = —- . ^ ^ -WTT’ (3) 

e + fi- e(p) - E(e) - Sp(e) 

where £(p) is the “bare” electronic dispersion, 
while the total self-energy completely neglects the 
interference between the Hubbard and additional 
interaction and is given by the additive sum of 
the local self-energy E(e) of DMFT and “exter¬ 
nal” self-energy Sp(e). This conserves the stan¬ 
dard structure of DMFT equations However, 
there are two important differences with standard 


DMFT. At each iteration of DMFT cycle we recal¬ 
culate the “external” self-energy Sp(£) using some 
approximate scheme for the description of “exter¬ 
nal” interaction and the local Green’s function is 
“dressed” by Ep(e) at each step of the standard 
DMFT procedure. 

In Fig. [T] we show the typical “skeleton” dia¬ 
grams for self-energy in DMFT-I-E. Here the first 
two terms are local DFMT self-energy diagrams 
due to Hubbard interaction, while two diagrams in 
the middle show contributions to self-energy from 
additional interaction (dashed interaction lines), 
and the last diagram (b) is a typical example of in¬ 
terference process which is neglected. Indeed, once 
we neglect such interference the total self-energy is 
defined as a simple sum of two contributions shown 
in Fig. UKa). 

As an effective Anderson impurity solver in our 
DMFT calculations we have always used the nu¬ 
merical renormalization group (T7j | , which allows to 
perform calculations at pretty low temperatures. 

For the self-energy due to disorder scattering 
produced by the Hamiltonian o we use below the 
simplest approximation neglecting the diagrams 
with “intersecting” interaction lines (like those in 
the fourth diagram of Fig. Ufa)), i.e. the so called 
self-consistent Born approximation, represented by 
the third diagram in Fig. [IJa). For the Gaussian 
distribution of site energies it is momentum inde¬ 
pendent and is given by: 

Sp(e)-)> ^G(£,p), (4) 

p 

where G{e, p) is the single-particle Green’s func¬ 
tion ©, while A is the strength of site energy dis¬ 
order. 

In the following we shall consider mainly the 
three-dimensional system with “bare” semi - el¬ 
liptic density of states (per unit cell and one spin 
projection), with the total bandwidth 2D, which 
is given by: 

No{e) = (5) 

In this case can directly demonstrate, that in 
DMFT-kE approximation disorder influence upon 
single - particle properties of disordered Hubbard 
model (both repulsive and attractive) is completely 
described by effects of general band widening by 
disorder scattering. Actually, in the system of 
self - consistent equations DMFT-I-E equations 
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[IlEITl both the “bare” band spectrum and 
disorder scattering enter only on the stage of cal¬ 
culations of the local Green’s function: 

G,, = ^G(e,p), (6) 

p 


where the full Green’s function G(£:,p) is deter¬ 
mined by Eq. while the self - energy due to 
disorder, in self - consistent Born approximation, 
is given by Eq. Then, the local Green’s func¬ 
tion takes the following form: 


fD 


G,, = 


de' - 


No{e') 


' —D £ ^ — s' — S(£) — /S?Gii 


-D 


Et-e' 


where we have introduced Et = e -I- /i — E(e) — 
A'^Gii- In the case of semi - elliptic density of 
states this integral can be calculated in analytic 
form, so that the local Green’s function reduces to: 


represents the density of states in the absence of 
interaction U widened by disorder. The density 
of states in the presence of disorder remains semi 
- elliptic, so that all effects of disorder scatter¬ 
ing on single - particle properties of disordered 
Hubbard model in DMFT-I-E approximation are 
reduced only to disorder widening of conduction 
band, i.e. to the replacement D —> He// ■ 

Within DMFT-I-S approach we can also inves¬ 
tigate the two-particle properties [l^, [2l|. After 
the general a naly sis, based on Ward identity de¬ 
rived in Ref. [^, we can show that the real part 
of dynamical (optical) conductivity in DMFT-I-E 
approximation is given by [l^ [2l| : 


Recr(a;) = ^ J de [/(£_) - /(e+)] x 
!]«(£+)- 


xRe 


1 - 


LO 

UJ 


Gii — 2 


Et - ^Ef - H2 


H2 


( 8 ) 


It can be easily seen that Eq. (jS]) represents one of 
the roots of quadratic equation: 


-1 

Gii — Et - 


( 9 ) 


reproducing the correct limit of Ga —>■ E^^ for 
infinitely narrow {D —>■ 0) band. Then we can 
write: 




G-d = £ + ^l- E{£) - A^Gu - —G,, = 

= £ + fj, — E(e)-(10) 


where we have introduced He// “ an effective half¬ 
width of the band (in the absence of electronic cor¬ 
relations, i.e. for U = 0) widened by disorder scat¬ 
tering: 


H 


eff 


= D\ 1-1-4 


H2' 


( 11 ) 


Now comparing ®, ® and m, we immediately 
see, that the local Green’s function can be written 
as: 


Gu = 


-Deff 


e + , 


Here 


^o(e) 


2 





( 12 ) 


(13) 


(14) 

where e is electronic charge, /(e±) — Fermi distri¬ 
bution with £± = £ ± ^ and 

^ORR(RA) q) _ ^ORRiRA) 

= hm-=-, 

g^-o g2 

(15) 

where the two-particle loops (see details in Ref. 
[Tsil ') q) contain all vertex correc¬ 

tions from disorder scattering, but does not in¬ 
clude any vertex corrections from Hubbard in¬ 
teraction. This considerably simplifies calcula¬ 
tions of optical conductivity within DMFT-I-E ap¬ 
proximation, as we have only to solve the single¬ 
particle problem determining the local self-energy 
E(e±) via the DMFT-l-E procedure, while non¬ 
trivial contributions from disorder scattering en¬ 
ter only via (w, q), which can be calcu¬ 

lated in some appropriate approximation, neglect¬ 
ing vertex corrections from Hubbard interaction. 
To be more specific, to obtain the loop contribu¬ 
tions (w, q), determined by disorder scat¬ 

tering, we can either use the usual “ladder” ap¬ 
proximation for the case of weak disorder, or fol¬ 
lowing Ref. [l^, we can use the direct general¬ 
ization of the self-consistent theory of localization 
(23 - [3nl| , which allows us to treat the case of strong 
enough disorder. In this approach conductivity 
is determined mainly by the generalized diffusion 
coefficient obtained from simple extension of self- 
consistency equation (28l®il| of this theory, which 
is to be solved in combination with DMFT-I-E pro¬ 
cedure [l8l| . 
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III. MOTT - ANDERSON TRANSITION 
IN DISORDERED SYSTEMS 

Below we present some of the most interesting 
results for repulsive Hubbard model at half-filling 
with semi - elliptic bare density of states (O with 
the bandwidth 2D [I^, which qualitatively is well 
suited to describe the three - dimensional case. 
Density of states below is given in units of number 
of states in energy interval for cubic unit cell of the 
volume a? (a is the lattice constant) and for one 
spin projection. Conductivity values are always 
given in natural units of /ha. 

A. Evolution of the density of states 

In the standard DMFT approximation density 
of states of repulsive Hubbard model at half-filling 
has a typical three-peak structure d, 0, [13 with 
pretty narrow quasiparticle (central) peak at the 
Fermi level and rather wide upper and lower Hub¬ 
bard bands situated at energies e ^ ±{7/2. As 
Hubbard repulsive interaction U grows quasipar¬ 
ticle band narrows within the metallic phase and 
disappears at Mott-Hubbard metal-insulator tran¬ 
sition at critical interaction value Uc2/‘2D « 1.5. 
At larger values of U we observe insulating gap at 
the Fermi level. 

In Fig. [5] we present our results [l3 for 
DMFT±S densities of states for typical strongly 
correlated metal with U = 2.5D, both in the ab- 
sense of disorder and for different values of disor¬ 
der scattering A, including strong enough values 
of disorder, which transforms correlated metal to 
correlated Anderson insulator. In metallic phase 
disorder scattering leads to a typical broadening 
and suppression of the density of states. 

Much more unusual is the the result obtained 
for U = 4.51?, typical for Mott insulator phase 
and shown on Fig. HDb). Here we observe the re¬ 
covery of the central peak (quasiparticle band) in 
the density of states with the increase of disorder, 
transforming Mott insulator to correlated metal or 
to correlated Anderson insulator. Similar density 
of states behavior for disordered Hubbard model 
was reported also in Ref. , using direct numer¬ 
ical DMFT calculations in finite lattices. 

Physical origin of such quite unexpected central 
peak restoration is evident. Controlling parameter 
of metal-insulator transition in DMFT is the ratio 
of Hubbard interaction U to bare bandwidth 2D. 
Introduction of disorder (in the absense of Hub¬ 
bard interaction) leads to new effective bandwidth 
2.Df,ff (cf. (fTT]) b growing with disorder. This leads 
to diminishing values of the ratio U/2Deff, which 
in its turn causes restoration of the quasiparticle 
band. 

More so, in complete accordance with analytic 




e/2D 


FIG. 2: Density of states of Hubbard model at half¬ 
filling for different disorder levels A [T^. (a) — cor¬ 
related metal with U = 2.5D. (b) — Mott insulator 
with U = 4.5D. Temperature T/2D=0.0005. 


arguments presented above, the density of states 
behavior in disordered Hubbard model with semi 
- elliptic actually demonstrates the universal de¬ 
pendence on disorder. This is clearly seen from 
Fig. 131 where we show properly normalized typical 
densities of states 2DeffN{e) in metallic (normal¬ 
ized interaction value {7/2i?e//=I.O) and insulat¬ 
ing phase (corresponding to U/2Deff=3.0) with¬ 
out disorder and for the typical value of disorder 
scattering A./2D =0.25. The densities of states 
in the absence and in the presence of disorder are 
actually described by the same (universal) depen¬ 
dences if expressed via properly normalized param¬ 
eters. 

In the absense of disorder one of the charac¬ 
teristic features of Mott-Hubbard metal-insulator 
transition is hysteresis behavior of the density of 
states, appearing with the decrease of U, start¬ 
ing from insulating phase 113 • Mott insula¬ 
tor phase remains (meta)stable down to rather 
small values of U deep within the correlated metal 
phase and metallic phase is restored only at about 
Uci/2.D K, 1. Corresponding interval of interac¬ 
tion parameter Uci < U < Uc 2 represents a co- 
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FIG. 3: Universal dependence of properly normalized 
density of states on normalized energy e/2De// in Hub¬ 
bard model for different disorder levels A. (a) — cor¬ 
related metal {U/2Def f=1.0) with no disorder and for 
Af2D =0.25. (b) — Mott insulator ([//2De//=3.0) 
without disorder and for Aj2D =0.25. Temperature 
r/2De//=0.0009. 


existence region of metallic and Mott insulating 
phases, where, from a thermodynamic point of 
view, metallic phase is more stable [y,[^[^. Such 
hysteresis in density of states behavior is observed 
also in the presence of disorder Blii. 


B. Optical conductivity: Mott-Hubbard and 
Anderson transitions 

In the absence of disorder our calculations repro¬ 
duce conventional DMFT results with op¬ 

tical conductivity characterized by typical Drude 
peak at low frequencies and wide maximum at 
about ui ^ U, which corresponds to optical transi¬ 
tions to the upper Hubbard band. As U grows 
Drude peak is suppressed and disappears com¬ 
pletely at Mott transition. Introduction of disor¬ 
der leads to qualitative changes of the frequency 
dependence of optical conductivity. 

Fig. |4l[a) shows the real part of optical conduc¬ 
tivity of Hubbard model at half-filling for differ¬ 




FIG. 4: Real part of optical conductivity of Hubbard 
model at half-filling for different disorder levels A [l^ . 
(a) - typical correlated metal with U = 2.5D. Curves 
1,2, - metallic phase, curve 3 corresponds to the mo¬ 
bility edge (Anderson transition), curves 4,5 - corre¬ 
lated Anderson insulator, (b) - typical Mott insula¬ 
tor with U = 4.5D. Curves 1,2 correspond to Mott 
insulator, curve 3 - mobility edge (Anderson transi¬ 
tion), curves 4,5 - correlated Anderson insulator. Inset 
shows magnified low frequency region. Temperature 
r/2D=0.0005. 


ent disorder levels A and U = 2.5D typical for 
correlated metal. Transitions to the upper Hub¬ 
bard bands at energies lo ^ U are almost unob¬ 
servable. However it is clearly visible that metallic 
Drude peak typically centered at zero frequency is 
broadened and suppressed by disorder, gradually 
transforming into a peak at finite frequency be¬ 
cause of Anderson localization effects. Anderson 
transition takes plase at Ac ~ Q.IAD (correspond¬ 
ing to the curve 3 on all figures here). Notice that 
this value explicitly depends on value of the cutoff 
in the equation for the generalized diffusion coef¬ 
ficient, which is defined upto the coefficient of the 
order of unity [2^, [2^. Naive expectations can lead 
to conclusion that narrow quasiparticle band at the 
Fermi level (formed in a strongly correlated metal) 
may be localized much more easily than the usual 
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conduction band. However, these expectations are 
wrong and the band localizes only at rather large 
disorder Ac ^ I? of the order of conduction band 
width ~ 2D. This is in qualitative agreement with 
the results for localization transition in two-band 
model 

In the DMFT-I-S approach critical disorder 
value Ac does not depend on U as interaction ef¬ 
fects enter here only through A'E^^{u!) —>■ 0 for 
w —^ 0 (for T = 0, e = 0), so that the influ¬ 
ence of interaction at w = 0 just disappears. In 
fact this is the main shortcoming of DMFT-I-E ap¬ 
proach originating from the neglect of interference 
effect between interaction and impurity scattering. 
Significant role of these interference effects is actu¬ 
ally well known for a long time However, 

the neglect of these effects allows us to perform 
the reasonable physical interpolation between two 
main limits - that of Anderson transition because 
of disorder and Mott-Hubbard transition because 
of strong correlations. 

On Fig. HKb) we show the real part of opti¬ 
cal conductivity of Mott-Hubbard insulator with 
U = 4.51? for different disorder levels A. In the in¬ 
set we show low frequency data, demonstarting dif¬ 
ferent types of conductivity behaviour, especially 
close to Anderson transition and within the Mott 
insulator phase. On the main part of the figure 
contribution to conductivity from transitions to 
upper Hubbard band at about a; ~ 17 is clearly 
seen. Disorder growth results in the appearance of 
finite conductivity for the frequencies inside Mott- 
Hubbard gap, correlating with the restoration of 
quasiparticle band in the density of states within 
the gap as shown in Fig. [2Kb). This conductivity 
for A < Ac is metallic (finite in the static limit 
w = 0), and for A > Ac at low frequencies we 
get Reg(a;) uj'^, which is typical for Anderson 
insulator [M. [28| - I^ . 

A bit unusual is the appearance in ReCT(a;) of 
a peak at finite frequencies even in the metallic 
phase. This happens because of importance of lo¬ 
calization effects. In the “ladder” approximation 
for $c^"^(a;,q) which neglects all localization ef¬ 
fects we obtain the usual Drude peak at w = 0 
[T^, while taking into account localization effects 
shifts the peak in ReCT(a;) to finite frequencies. 

Above we presented the data for conductivity 
data obtained for the case of increase of U from 
metallic to Mott insulator phase. As U decreases 
from Mott insulator phase we observe hysteresis 
of conductivity in coexistence region defined (in 
the absense of disorder) by inequality Uci < U < 
Uc 2 - Hysteresis of conductivity is also observed in 
the coexistence region in the presence of disorder. 
More details can be found in Refs. El, El. 

In general, the picture of conductivity behav¬ 
ior obtained in DMFT-I-S approximation is rather 
rich, demonstarting both Mott - Hubbard tran- 



FIG. 5: Phase diagram of disordered Hubbard model 
El. Continuous curves are Mott insulator phase 
boundaries 17ci,c2(A) obtained from analytical esti¬ 
mate of Eq. G3, different symbols represent results for 
these boudaries obtained from calculations from den¬ 
sity of states and optical conductivity. Line of Ander¬ 
son transition is given by Ac = 0.37. 

sition due to strong correlations and disorder in¬ 
duced Anderson (localization) transition. The 
complicate behavior under disordering is essen¬ 
tially determined by two - particle Green’s func¬ 
tion behavior and does not shown a kind of univer¬ 
sality, demonstrated above for the single - particle 
density of states. 


C. Phase diagram of disordered Hubbard 
model at half-filling 

Phase diagram of repulsive disordered Hubbard 
model at half-filling was studied in Ref. [Slj , using 
direct DMFT numerics for lattices with finite num¬ 
ber of sites with random realizations of energies Ci 
in o, with subsequent averaging over many lat¬ 
tice realizations to obtain the averaged density of 
states and geometric mean local density of states, 
which allows to determine the critical disorder for 
Anderson transition. Below we present our results 
on disordered Hubbard model phase diagram ob¬ 
tained from density of states and optical conduc¬ 
tivity calculations in DMFT-t-S approach E^ - 
Calculated disorder - correlation strength 
(A, U) phase diagram at zero temperature is shown 
in Fig. |S] (actually calculations were performed 
at very low T/2D—0.0005). Anderson transition 
line Ac ~ 0.747? is defined as disorder strength for 
which static conductivity becomes zero at T = 0. 
Mott-Hubbard transition can be detected either 
from central peak disappearance in the density of 
states or from optical conductivity by observation 
of gap closing in the insulating phase or from static 
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conductivity disappearance in the metallic phase. 

We have already noticed that DMFT+S approx¬ 
imation gives universal {U independent) value of 
critical disorder Ac because of neglect of interfer¬ 
ence between disorder scattering and Hubbard in¬ 
teraction. This leads to the difference between the 
phase diagram of Fig. [5] and the one obtained by 
numerical simulations in Ref. [^ . At the same 
time the qualitative form of our phase diagram is 
highly nontrivial and qualitatively coincide with 
results of Ref. [3l| . Main difference is conservation 
of Hubbard bands in our results even in the limit 
of high enough disorder, while in the Ref. [Slj they 
just disappear. Phase coexistence region in Fig. 0 
slowly widens with disorder growth instead of van¬ 
ishing at some “critical” point as on phase diagram 
of Ref. [3l| . Coexistence boundaries (Mott insula¬ 
tor phase boundaries), obtained with decrease or 
increase of U, represented by curves C/ci(A) and 
t/c 2 (A) on Fig. [5l can actually be obtained from 
the simple equation: 

Ucl ,c2 (A) C/ci .c2 

~~n -= ~~n~’ 

^eff O 

where effective bandwidth in the presence of dis¬ 
order is calculated for C/ = 0 within self-consistent 
Born approximation o, m- Thus the bound¬ 
aries of coexistence region (which define also the 
boundaries of Mott insulator phase) are given by: 

I A^ 

Hc1,c 2(A) = C/cl,c2 Y 1 + 4-^ (17) 

which are shown in Fig. [S] by dotted and solid 
lines. Phase transition points detected from dis¬ 
appearance of quasiparticle peak as well as points 
following from qualitative changes of conductivity 
behaviour are shown in Fig. [S] by different sym¬ 
bols. These symbols demonstrate very good agree¬ 
ment with analytical results confirming the choice 
of the ratio (ITSl) as a controlling parameter of Mott 
transition in the presence of disorder. Thus, this 
transition is essentially controlled by simple band 
- widening effects due to disorder scattering, sim¬ 
ilarly to density of states behavior demonstrated 
above. 

Note that the values of normalized density of 
states 2Def fN{e) are universal along each of these 
boundaries, as well as along any curve in {A,U) - 
plane, determined by the equation: 

C/(A) = C/(0)y/i + 4^ (18) 

in accordance with our discussion on the universal 
dependence of the densities of states on disorder 
presented above. 

Essentially similar results were obtained for the 
density of states behavior, dynamic conductivity 


and phase diagram [T^ in the case of the conduc¬ 
tion band with “flat” density of states in the ab¬ 
sence of disorder and interactions, which qualita¬ 
tively corresponds to the two - dimensional case. 
This is not surprising, as both large enough dis¬ 
order and interactions transform the “flat” band 
into a kind of smeared semi - elliptic band. Some 
explicit examples of this kind of behavior will be 
presented below for the case of attractive Hubbard 
model. 


IV. ATTRACTIVE HUBBARD MODEL 
WITH DISORDER 

The studies of superconductivity in BCS - BEC 
crossover region attracts theorists for rather long 
time 0 and most important advance here was 
made by Nozieres and Schmitt-Rink [llj |. who 
proposed an effective approach to describe Tc 
crossover. Attractive Hubbard model is probably 
the simplest model allowing theoretical studies of 
BCS-BEC crossover m- This model was studied 
within DMFT in a number of recent papers [ 13 - 
[iojl . However only few results were obtained for the 
normal (non-superconducting) phase of this model, 
especially in disordered case. Similarly, there were 
practically no studies of two - particle properties, 
such as optical conductivity. Below we present a 
summary of our results obtained within DMET-I-E 
approach and make comparison with similar re¬ 
sults for repulsive Hubbard model. 


A. Density of states and optical conductivity 

In the special case of half - filled band (n = 1) 
the densities of states of attractive and repulsive 
Hubbard models just coincide (due to exact map¬ 
ping of these models onto each other). Thus, be¬ 
low we discuss the more typical case of quarter - 
filled band {n = 0.5). In Eig. [5] we show densi¬ 
ties of states obtained for T /2Z) = 0.05 for differ¬ 
ent values of attractive interaction {U < 0). Fig. 
IHKa) should be compared with Fig. EKb), where we 
present similar results for repulsive {U > 0) case. 
We can see that the densities of states close to the 
Fermi level drop with the growth of 17, both for 
attraction (Fig. UKa)) and repulsion (Fig. HKb)), 
but significant growth of |t/| in repulsive case leads 
only to vanishing quasiparticle peak, so that the 
density of states at the Fermi level becomes prac¬ 
tically independent of U, while in attractive case 
the growth of |C/| leads to superconducting pseu¬ 
dogap opening at the Fermi level (curves 5, 6 in 
Fig. [HKa)) and for \U\/2D > 1.2 we observe the 
full gap opening at the Fermi level (curves 7-9 in 
Fig. [IK a)). This gap is not directly related to the 
emergence of superconducting state, but is due to 






0,5 



e/2D 



FIG. 6: Densities of states and for different values of 
Hubbard attraction (a) and repulsion (b). Tempera¬ 
ture T/2D = 0.05. 









FIG. 7: Optical conductivity for different values of at¬ 
tractive Hubbard attraction. Temperature T/2D = 
0.05. 


the appearance of preformed Cooper pairs at the 
temperatures larger than superconducting transi¬ 
tion temperature (which is lower, than the tem¬ 
perature T/2D = 0.05 used in our calculations). 
Here we actually observe the important difference 
between attractive and repulsive cases — in case 
of repulsion deviation from half - filling leads to 
metallic state for arbitrary values of U and insu¬ 
lating gap at large U opens not at the Fermi level. 

This picture of density of states evolution with 
the growth of |C/| is also supported by the behavior 
of optical conductivity shown in Fig. [T] We see 
that the growth of |C/| leads to the replacement of 
Drude peak at zero frequency (curves 1-3 in Fig. 
[7]) by pseudogap dip (curves 5, 6 in Fig. [7]) and 
wide maximum of conductivity at finite frequency, 
connected with transitions across the pseudogap. 
The further increase of |?7| leads to opening of the 
full gap in optical conductivity due to formation of 
Cooper pairs (curves 7-9 in Fig. Q. 

In Fig. [5] we present the evolution of the density 
of states and optical conductivity with changing 
disorder. At weak enough attraction {\U\/2D = 


FIG. 8: Evolution of the density of states (left pan¬ 
els) and optical conductivity (right panels) with disor¬ 
der for different values of U {\U\/2D = 0.8 - (a),(b); 
|[/|/2D = 1 - (c),(d); |[/|/2D = 1.6 - (e),(f)). 


0.8, Fig. l8Da),(b)), the growth of disorder just 
widens the density of states. Disorder effectively 
masks peculiarities of the density of states due 
to correlation effects. In particular, quasiparticle 
peak and the “wings” due to upper and lower Hub¬ 
bard bands present in Fig. IS^a) in the absence of 
disorder completely vanish at strong enough dis¬ 
order. Evolution of optical conductivity with the 
growth of disorder A, shown in Fig. HKb), is in 
general agreement with the evolution of density of 
states. Weak enough disordering (curves 1-3 in 
Fig. IHKb)), leads to some growth of static conduc¬ 
tivity, which is connected with suppression of cor¬ 
relation effects at the Fermi level (curves 1-3 in Fig. 
Ha). The further growth of disorder leads to sig¬ 
nificant widening of the band and the drop of the 
density of states (curve 5 in Fig. |lKa),(b)), which 
leads to the drop of static conductivity. Finally, 
the growth of disorder leads to Anderson localiza¬ 
tion which takes place at A/2Z1 = 0.37 for T = 0 
[ 3 . However, here we consider the case of high 
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enough temperature T/2D = 0.05, so that static 
conductivity (see curves 6, 7 in Fig. |HKb)) always 
remains finite, though the localization behavior is 
also clearly seen and cr(a;) ^ At larger value 
of attractive interaction \U\/2D = 1 the evolution 
of the density of states and optical conductivity 
is more or less similar (Fig. [8Kc),(d) ). However, 
in the absence of disorder we observe here Cooper 
pairing pseudogap in the density of states, while 
disorder leads to its suppression, leading both to 
the growth of the density of states at the Fermi 
level and related growth of static conductivity. Fi¬ 
nally, at still larger attraction \U\/2D = 1.6 (Fig. 
|SKe),(f)) in the absence of disorder there is the real 
Cooper pairing gap in the density of states. This 
gap is also evident in optical conductivity. With 
the growth of disorder Cooper pairing gap both 
in the density of states and conductivity becomes 
narrower (curves 2-5). Further growth of disor¬ 
der leads to complete suppression of this gap and 
restoration of metallic state with finite density of 
states at the Fermi level and finite static conduc¬ 
tivity. This closure of Cooper gap is obviously 
related to the effective growth of the conduction 
bandwidth 2Zle//, which leads to the lowering of 
\U\/2Deff ratio, which actually controls the for¬ 
mation of Cooper gap. Situation here is similar 
to the closure of Mott gap by disorder in repulsive 
Hubbard model discussed above [l^. However, at 
large disorder (curve 7 in Fig. HKf)) we clearly ob¬ 
serve localization behavior, so that the growth of 
disorder at T = 0 will first lead to metallic state 
(the closure of Cooper pairing gap), while the fur¬ 
ther growth of disorder will induce Anderson metal 
- insulator transition. Similar picture is observed 
for large positive U at half-filling (n = 1) (T^ . 
where the growth of disorder leads to Mott insula¬ 
tor - correlated metal - Anderson insulator tran¬ 
sition. 

Let us now demonstrate the universality of dis¬ 
order dependence of the density of states as an 
example of the most important single - particle 
property. Let us concentrate on the most typical 
case of the density of states evolution shown in 
Fig. |S] (a). We can easily convince ourselves, that 
this evolution is only due to the general widening 
of the band due to disorder (cf. (ITT]) 1. as all the 
data for the density of states fit the same universal 
curve replotted in appropriate new variables, with 
all energies (and temperature) normalized by the 
effective bandwidth by replacing D —>• 7?e//, as 
shown in Fig. IHKa), in complete accordance with 
results obtained above in repulsive Hubbard model 
for semi - ellipric band. 

In the case of initial (“bare”) conduction band 
with flat density of states, there is no complete uni¬ 
versality, as is seen from Fig. Elb) for low enough 
values of disorder. However, for large enough dis¬ 
orders the dashed curve in Fig. EKb) practically 



FIG. 9: Universal dependence of the density of states 
on disorder in attractive Hubbard model: (a) — the 
model of semi - elliptic “bare” density of states, (b) 
— the model of flat “bare” density of states. 

coincides with universal curve for the density of 
states shown in Fig. [9Ka). This reflects the simple 
fact, that at large enough disorders the flat den¬ 
sity of states is effectively transformed into semi - 
elliptic one [13 . 

B. Generalized Anderson theorem 

Superconducting transition temperature in 
general is not a single — particle characteristic 
of the system. Cooper instability, determining 
Tc is related to divergence of two - particle loop 
in Cooper channel. In the weak coupling limit, 
when superconductivity is due to the appearance of 
Cooper pairs at disorder only slightly influences 
superconductivity with s-wave pairing [4l|, 1 ^ . 
This is the essence of the so called Anderson theo¬ 
rem and changes of are due only to the relatively 
small changes of the density of states at the Fermi 
level induced by disorder. 

In region of BCS - BEC crossover and in the 
strong coup ling region Nozieres - Schmitt-Rink ap¬ 
proach [ll| assumes, that corrections due to strong 
pairing attraction significantly change the chemical 
potential of the system, while possible correction 
due to this interaction to Cooper instability con¬ 
dition can be neglected, so that we can always use 
the weak coupling (ladder) approximation. Then 
the condition of Cooper instability in disordered 
Hubbard model takes the form: 

l = -|t^lxo(g = 0,cc™ = 0) (19) 

where 

Xo(g = 0,c.™ = 0) = -T^^$pp.(£„) (20) 

n pp' 

represents the two - particle loop (susceptibility) 
in Cooper channel “dressed” only by disorder scat¬ 
tering, and $pp'(en) is the averaged two - particle 
Green’s function in Cooper channel {tOm = 27rmT 









and En = 7rr(2n + 1) are the usual Boson and 
Fermion Matsubara frequencies). 

Using the exact Ward identity, derived in Ref. 
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G(e„, p) - G(-e„, -p) = 


P' 


( 21 ) 


where G(e„,p) is the impurity averaged (but not 
containing Hubbard interaction corrections!) sin¬ 
gle - particle Green’s function, we can show (2^ 
that Cooper susceptibility (l20ll is given by: 


Xo{q = 0,0;^ = 0) = 

_ Ep P) - Sp G{-en, p) 

n 

Epg(g»,p) 


( 22 ) 


After the standard summation over Matsubara fre¬ 
quencies [15 we get: 


Xo{q = 0, Wm = 0) = 


1 

47ri 


7-00 e 22 



de 


Noje) 

2s 



(23) 


where No{e) is the density of states {U = 0) renor¬ 
malized by disorder scattering. In Eq. (1^ the 
energy e origin is at the chemical potential. If the 
origin of energy is shifted to the middle of conduc¬ 
tion band we have to replace e —>■ e — /i, and the 
condition of Cooper instability (ITOl) leads to the 
following equation for Tc'. 




deNo{s) 


£ — 


(24) 


The chemical potential of the system at differ¬ 
ent values of U and A now should be determined 
from DMFT-I-E calculations, i.e. from the stan¬ 
dard equation for the number of electrons (band¬ 
filling), determined by Green’s function given by 
Eq. ([3]), which allows us to find Tc for the wide 
range of model parameters, including the BCS- 
BEC crossover and strong coupling regions, as well 
as for different levels of disorder. This is the gist 
of Nozieres - Schmitt-Rink approximation — in 
the weak coupling region superconducting transi¬ 
tion temperature is controlled by the equation for 
Cooper instability (El, while in the strong cou¬ 
pling limit it is determined by the temperature of 
Bose - Einstein condensation, which is controlled 
by chemical potential. Then the joint solution of 



FIG. 10: Dependence of superconducting critical tem¬ 
perature on attractive interaction strength. Black 
squares, white circles and white squares show the re¬ 
sults of Refs. [13],[15,[15 respectively for quarter- 
filled band with n — 0.5. Stars represent the results 
obtained numerically from the criterion of instability 
of the normal phase in Ref. [^ . Filled circles show Tc 
in Nozieres - Schmitt-Rink approximations, combined 
with DMFT [^. Continuous black curve represents 
the result of BCS theory. 


Eq. (1331) and equation for the chemical potential 
guarantees the correct interpolation for Tc through 
the region of BCS-BEC crossover. In the absence 
of disorder this combination of Nozieres - Schmtt- 
Rink approximation with DMFT produces the re¬ 
sults for the critical temperature, which, as shown 
in Fig. [TUI are almost quantitatively close to exact 
results, obtained W direct numerical DMFT cal¬ 
culations [ 13 , [ 13 , Hi, nil , but demands much less 
numerical efforts. 

Eq. (lUTl) demonstrates, that Cooper instability 
depends on disorder only through the disorder de¬ 
pendence of the density of states iVo(e), which is 
the main statement of Anderson theorem. Within 
Nozieres - Schmitt-Rink approach Eq. dm) is con¬ 
served also in the region of strong coupling, when 
the critical temperature is determined by BEC 
condition for compact Cooper pairs. However, the 
chemical potential fi, entering Eq. dm, may signif¬ 
icantly depend on disorder. In DMFT-I-E approx¬ 
imation this dependence of chemical potential (as 
well as any other single - particle characteristic) in 
the model with semi - elliptic density of states is 
only due to disorder widening of conduction band. 
In this sence both in BCS - BEC crossover region 
and in the strong coupling limit a kind of general¬ 
ized Anderson theorem actually holds and Eq. (IMl) 
leads to universal dependence of Tc on disorder, 
due to the change of D —>• Def /. Such universality 
is fully confirmed by direct numerical calculations 
of Tc in this model, performed in Ref. [^ . 

In Fig. [TT]we present the dependence of Tc (nor¬ 
malized by the critical temperature in the absence 
of disorder Tco = Tc{A — 0)) on disorder for differ- 
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FIG. 11: Dependence of superconducting transition 
temperature on disorder for different values of Hub¬ 
bard attraction U: (a) — semi - elliptic band, (b) — 
fiat band. 



FIG. 12: Universal dependence of superconducting 
critical temperature on Hubbard attraction U for dif¬ 
ferent disorder levels: (a) — semi - elliptic band. 
Dashed curve represent BGS dependence in the ab¬ 
sence of disorder, (b) — flat band. Dashed line rep¬ 
resents similar dependence for semi - elliptic band for 
A = 0. 


ent values of pairing interaction U for both mod¬ 
els of initial semi - elliptic density of states (Fig. 
[TlT aH and for the case of flat density of states (Fig. 
inrb)). Qualitatively the evolution of Tc with dis¬ 
order is the same for both models. In the weak 
coupling limit {U/2D <C 1) disorder slightly sup¬ 
presses Tc (curves 1). At intermediate couplings 
{U/2D ^ 1) weak disorder increases Tc, while the 
further growth of disorder suppresses the critical 
temperature (curves 3). In the strong coupling 
region {U/2D 1) the growth of disorder leads 

to significant increase of the critical temperature 
(curves 4,5). However, this rather complicate de¬ 
pendence of Tc on disorder is actually completely 
determined simply by disorder widening of the ini¬ 
tial {U = 0) conduction band, demonstrating the 
validity of the generalized Anderson theorem for all 
values of U. In Fig. [12] curve with octagons show 
the dependence of the critical temperature Tc/2D 
on coupling strength U/2D in the absence of disor¬ 
der (A = 0) for both models of initial conduction 
bands (semi - elliptic — Fig. fT^ al and flat — 
Fig. HUb)). In both models in the weak coupling 


region superconducting transition temperature is 
well described by BCS model (in Fig. (UKa) the 
dashed curve represents the result of the solution 
of BCS model, with Tc determined by Eq. (l24ll . 
with chemical potential independent of U and de¬ 
termined by quarter - filling of the “bare” band), 
while in the strong coupling region the critical tem¬ 
perature is determined by Bose - Einstein conde- 
sation of Cooper pairs and drops as /U with the 
growth of U (inversely proportional to the effective 
mass of the pair), passing through the maximum - 
at U/2Dcff ^ 1. The other symbols in Eig. IHKa) 
show the results for Tc obtained by combination of 
DMFT-l-I] and Nozieres - Schmitt-Rink approxi¬ 
mations for the case of semi - elliptic band. We can 
see, that all data (expressed in normalized units of 
U/2Dcff and Tc/2Dcff) ideally fit the universal 
curve, obtained in the absence of disorder. For the 
case of flat band, results of our calculations are 
shown in Fig. [T^ b'l and we do not observe the 
complete universality — data points, correspond¬ 
ing to different degrees of disorder slightly deviate 
from the curve, obtained in the absence of disorder. 
However, with the growth of disorder the flat den¬ 
sity of states gradually transforms to semi - ellip¬ 
tic and our data points move towards the universal 
curve, obtained for semi - elliptic case and shown 
by the dashed curve in Fig. UlTb), confirming the 
validity of the generalized Anderson theorem also 
in this case. 


C. Ginzburg — Landau coefficients 

Universal dependence on disorder is also ob¬ 
served for the coefficients of Ginzburg - Landau 
expansion A (homogeneous quadratic term of the 
expansion) and B (fourth-order term), related to 
Cooper - channel vertices with the sum of incom¬ 
ing (outgoing) momenta q = Q. Coefficient A is 
given by 

A(T) = Xo(q = 0, T) -Xo(q = 0, TJ, (25) 

where xo(9 = 0,T) is Cooper susceptibility (1^ . 
and subtraction of Xo(q = 0,Tc) guarantees the 
zero value of A(T = Tc). Using (fTOl) to determine 
Xo(q = 0,Tc) and ([121) for Xo(<? = 0,r), we get: 

SO that the coefficient A{T) reduces to zero for T —?> 
Tc, and is written as: 

A{T) = a(T-Tc). (27) 

For the case of “bare” band with semi - elliptic 
density of states the dependence of a on disorder 
is related only to the general widening of the band 
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coefficients A and B were otained using the exact 
Ward identities and remain valid also in the limit of 
strong disorder (Anderson localized phase), when 
both. A and B depend on disorder also only via 
the effective bandwidth D^ff. 

This universal dependence on disorder (due only 
to the replacement D —^ A*e//) is reflected also 
in the specific heat discontinuity at the transitio 
temperature, which is determined by coefficients a 
and B : 


Cs{T,)-Cn{Tc) = Tj^ (29) 

To determine the coefficient C in gradient term 
of Ginzburg - Landau expansion we need the 
knowledge of nontrivial of g-dependence of Cooper 
vertex [26j |. which is essentially changed by disor¬ 
der scattering. In particular, the behavior of coef¬ 
ficient C is qualitatively changed at Anderson lo¬ 
calization transition [^. Thus, the coefficient C 
is basically determined by two - particle chara- 
teristics of the system and does not demopnstrate 
the universal dependence on disorder due only to 
changes of the effective bandwidth. 


D. Number of local pairs 


FIG. 13: Universal dependence of Ginzburg - Landau 
coefficients a (a) and B (b) on Hubbard attraction for 
different disorder levels. Dotted line with black squares 
shows the case of A = 0. 


by disorder, i.e. is completely described by the 
replacement D —^ D^ff. Thus, in the presence of 
disorder we obtain the universal dependence of a 
on U (normalized by D^ff), shown in Fig. fT^ . 

Ginzburg - Landau coefficient B is determined 
by the “loop” diagramm with four Gooper vertices 
[^ . After rather complicated analysis, to be pre¬ 
sented elsewhere, which is based on some gener¬ 
alizations of Ward identity (EU, it can be shown 
exactly, that B is given by: 


de 

7-00 4(e-/i)3 



S — 

2T 


{e-fi)/2T \ 

I 


Noie) 


(28) 

Thus, the dependence of coefficient B on disorder, 
similarly to A, is determined only by the density 
of states No (e) renormalized (widened) by disorder 
and the chemical potential /i. Then, in the case of 
semi - elliptic density of states the dependence of 
B on disorder is reduced to the simple replacement 
D —^ Deff, so that in the presence of disorder we 
obtain again the universal dependence of B on t/, 
shown in Fig. 113b . 

It should be noted that Eqs. (E51) and (E51) for 


Disorder in attractive Hubbard model also leads 
to the suppression of the number of local pairs 
(doubly occupied sites). The average number of 
local pairs is determined by local (single site) pair 
correlation function < n-j-nj, >, which in the ab¬ 
sence of disorder grows with the increase of Hub¬ 
bard attraction U from < n-i-nj, >=< n-j- >< 
rif >= n^/4 for U/2Deff <C 1 to < n-i-nj, >= n/2 
for U/2Deff ^ 1, when all electrons become 
paired. In our calculations n=0.5 (quarter - filled 
band), so that n/2=0.25, while n^/4=0.0625. The 
growth of Deff with disorder leads to an effective 
suppression of the parameter U/2Deff and cor¬ 
responding suppression of the number of doubly 
occupied sites. In Fig. fldl al we show the dis¬ 
order dependence of the number of doubly occu¬ 
pied sites for three different values of Hubbard at¬ 
traction. In all cases the growth of disorder sup¬ 
presses the number of doubly occupied sites (local 
pairs). Similarly to Tc, the change of the num¬ 
ber of local pairs with disorder can be attributed 
only to the change of the effective bandwidth (HU 
with the growth of disorder. In Fig. [HKb) the 
curve with black squares shows the dependence 
of the number of doubly occupied sites on attrac¬ 
tive interaction in the absence of disorder at tem¬ 
perature T/2D = 0.0586. This curve is actu¬ 
ally universal — the dependence of the number 
of local pairs < n-j-nj, > on the scaled parame¬ 
ter U/2Deff with appropriately scaled tempera¬ 
ture T/2Deff = 0.0586 in the presence of disorder 
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FIG. 14: Dependence of the number of local pairs on 
disorder for different values of Hubbard attraction (a) 
and universal dependence on disorder, expressed via 
normalized variables at fixed value of U/2D = 1 (b). 

is given by the same curve, which is shown by by 
circles representing data obtained for five different 
disorder levels, shown in Fig. [T^ bl for the case of 
U/2D = 1. 

V. CONCLUSION 

In this paper, in the framework of DMFT+E 
generalization of dynamical mean field theory [13: 
we have studied and compared disorder effects in 
both repulsive and attractive Hubbard models. We 
examined both the problem of Mott - Hubbard 
and Anderson metal-insulator transitions in repul¬ 
sive case, and BCS-BEC crossover region of attrac¬ 
tive Hubbard model. We also performed extensive 
calculations of the densities of states and dynamic 
(optical) coductivity for the wide range of inerac- 
tions U and at different disorder levels A, demon¬ 
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strating similarities and dissimilarities between re¬ 
pulsive and attractive cases. 

We have shown analytically for the case of con¬ 
duction band with semi - elliptic density of states 
(which is a good approximation for three - dimen¬ 
sional case) in DMFT-I-E approximation disorder 
influences all single - particle properties (e.g. den¬ 
sity of states) in a universal way — all changes of 
these properties are due only to disorder widening 
of the conduction band. In the model of conduc¬ 
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Naturally, no universal dependences on disor¬ 
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like optical conductivity, where vertex corrections 
due to disorder scattering become very important, 
leading to new physics, like that of Anderson tran¬ 
sition. 

Overall, the use of DMET-I-E approximation to 
analyze the disorder effects in Hubbard model was 
shown to produce reasonable results for the phase 
diagram for repulsive case, as compared to exact 
numerical simulations of disorder in DMFT, den¬ 
sity of states behavior and optical conductivity in 
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role of approximations made in DMET-I-E, such as 
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and correlation effects, deserves further studies. 
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